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Depression Disorder Integrated Model (DDIM) 

A method of choosing treatment based on selection (after outcome predictions 
are made) and then a method for generating a profile for patient monitoring 
patient response to specific treatment based on a method for prediction 
treatment specific response patterns. 



Diagnosis and other baseline data are collected and transmitted as input to 
the "Outcome Predictor." The outcome predictor has been trained on other 
(prior) patients baseline and outcome data. Predictions and confidence limits 
are generated for a number of treatment options. 

Treatment selection suggestion is made by the system for evaluation and 
consideration by a clinician, who has the final responsibility for the 
treatment selection. 

The baseline data and the chosen treatment are-input to the Pattern 
Predictor 
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Training Cycle for Each Treatment Group 
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Figure 3-2: Training cycle for each treatment group. Each of the two models with separate sets 
of parameters, but with a same architecture, was trained by cycling through individual data within 
the respective treatment group. After each cycle, the cost function (which reflects the degree of fit 
of the model predictions to the actual data) was evaluated to determine the completion of training. 
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Figure 3*4: The annotated second order differential equation used to model the pattern of 
recovery. 
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Figure 3-5: Direct effects of treatment are either immediate (step function) or delayed (sig- 
moid function). Delays are estimated by treatment from the patient data, using an optimization 
procedure. 
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Figure 3-6: Direct effects and interactions in the recovery model, u; represents the strength 
of the immediate effect of treatment on symptom node i; m represents the strength of the delayed 
effect of treatment on symptom i; and wji and w {j represent the interaction between the symptoms: 
the strength of the effect of symptom i on symptomj and the strength of the effect of symptom j 
on symptomi, respectively. 
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Figure 3-8: A schematic description of the cost function L. The formula inside 
the integration has two terms. By minimizing the first term, discrepancies between 
estimated and actual patterns of recovery are minimized. The second term is a 
constraint term which states that the differential equation must hold. 
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(a) CBT patient 180 
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(b) CBT patient 195 
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(c) DMI patient 155 
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Figure 3-10: Individual patterns of recovery. Plots, except the one at the bottom right, 

show patterns of recovery. Solid lines show actual patterns, dotted lines show predicted patterns. 
Numbers shown at vertical axis are scaled such that the possible maximum symptom factor value 
yields 1.0. The plot at the bottom right shows the error (L) on the ordinate axis plotted against 
number of training cycles on the abscissa. 
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(a) Predicted patterns of recovery by shunting equations 
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(b) Predicted patterns of recovery by second order equations 
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(a) Predicted and actual mean patterns of recovery (CBT) 
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(b) Predicted and actual mean patterns of recovery (DMI) 

Figure 3-11: Predicted and actual recovery patterns of patient data mean, (a) mean of six 
CBT responders (b) mean of six DMI responders. In both (a) and (b), all plots except the sub- 
plots at the bottom right, show mean patterns of recovery. Solid lines show actual mean patterns, 
dotted lines show predicted mean patterns. Numbers shown at vertical axis are scaled such that 
the maximum possible symptom factor value yields 1.0. The plot at the bottom right in both (a) 
and (b) shows the error L plotted on the ordinate axis against number of training cycles on the 
abscissa. Pat 0 indicates the patient averaged patient data were used. 
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Figure 3-13: Predicted CBT and DMI temporal response sequence of symptom improvement. 
The vertical axis shows the mean half reduction time in weeks; the horizontal axis has no meaning. 
Symptom names are placed vertically at their mean half reduction time. Significant difference 
(p < 0.05 after rounding) between half reduction times of energy and early sleep disturbance. There 
is a trend (p < 0.10 after rounding) for a split between cognitions and work in CBT responders. In 
DMI there were no significant differences (or trends) in the sequence. 
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Figure 31o: Graphic representation of the sequence of symptom factors in recovery with 
Cognitive Behavioral Therapy treatment for the second order system. Vertical positions of the 
symptoms represent half-way-reduction time, arrows represent strong impacts and interactions 
and corresponding numbers indicate the strength of the impact or interaction. See text for formula 
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Figure 3*16: Graphic representation of the sequence of symptom factors in recovery with 
Desipramine treatment for the second order model. Vertical positions of the symptoms represent 
half-way-reduction time, arrows represent strong impacts and interactions, and corresponding num- 
bers indicate the strength of the impact or interaction. Dotted arrows show the interactions that 
operate in loops. See text for formula. 
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Table 4.1: Prior attempts to predict outcome from Hamilton Depression Rating Scale scores 
(symptoms) and severity (total HDRS) of depressive symptoms. Independent variables are the 
HDRS items described under the column Symptom. Numbers are indices into Table 4.2 where the 
citations are listed and indicate that the symptom at baseline was not found to be a significant pre- 
dictor of outcome. Underlined numbers are also indices into Table 4.2, but indicate the symptom at 
baseline was found to predict outcome. Blank entries indicate that the significance of the symptom 
to predict outcome was not reported. (+) indicates greater intensity predicted better response; (-) 
indicates greater intensity predicted poorer response. See Table 4.2 for key to citation references. 
Other=S-adenosyl methionine (Carney et al., 1986), ECT or unspecified antidepressant medica- 
tion (Hinrichsen k Hernandez, 1993), unspecified antidepressant medication (Katon et al., 1994); 
Tricyclic=tricyclic antidepressant medications; MAOI=monoamine oxidase inhibitor antidepres- 
sant medications; SSRI=selective serotonin reuptake inhibitors; ECT=electroconvulsive therapy; 
Ami=amitriptyline; Nor=nortriptyline; IMI=imipramine; DMI=desipramine; Clo=clomipramine; 
Lev=levoprotiline; Map=maprotiline; Phe=phenelzine; Tran=tranylcypromine; Flu=fluoxetine; 
lPT=interpersonal therapy; CBT=cognitive behavioral therapy. 
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Figure 4-1: Nonlinear mapping of backpropagation. Each hidden node finds a direction in the 
input space (illustrated by an arrow perpendicular to a small square piece) to which the output 
is sensitive to. The output of each hidden node goes through a nonlinear output function before 
being weighted and summed at the output node. 
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Figure 4*2: Schematic representation of the effect of normalizing transformations on reducing 
nonlinearity of score-to-output relationships (or skewedness of distributions). In the transformation, 
the area under the curve is preserved. The transformation redistributes the position of the data 
values along the x-axis in order to preserve the areas under the curve between adjacent scores values 
while redistributing these data to best approximate a normal distribution. Equal areas under the 
curve between percentiles map to equal areas under the curve in the new distribution. 
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